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Abstract 

In order to depict the deviation of quantum time evolution in open systems from Markovian 
processes, different measures have been presented. We demonstrate that the measure proposed by 
Breuer, Laine and Piilo [Phys. Rev. Lett. 103, 210401 (2009)] and the two measures proposed 
by Rivas, Huelga and Plenio [Phys. Rev. Lett. 105, 050403 (2010)] have exactly the same non- 
Markovian time-evolution intervals and thus are really equivalent each other when they apply to 
open two-level systems coupled to environments via Jaynes-Cummings or dephasing models. This 
equivalence implies that the three measures in different ways capture the intrinsical characters 
of non-Markovianty of quantum evolutional processes. We also show that the maximization in 
the definition of the first measure can be actually removed for the considered models without 
influencing the sensibility of the measure to detect non-Markovianty. 
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I. INTRODUCTION 



The evolution of open quantum systems can be divided into two basic types, Markovian 
and non- Markovian processes For Markovian processes, the correlation time between 
the system and environment is considered to be infinitesimally small so that the dynamical 
map does not carry any memory effects, leading to a monotonic flow of the information from 
system to environment. In contrast, non-Markovian processes with memory have different 
dynamical traits which give rise to the back flow of information from environment to the 
system 

as 

. Recently, people found that non-Markovian processes can lead to distinctly 
differ effects on decoherence and disentanglement Q Q of open systems compared to 
Markovian processes. Many relevant physical systems, such as quantum optical system 

s, 



quantum dot 



6|, color-core spin in semiconductor [7|, could not be described simply by 



Markovian dynamics. In addition, quantum chemistry [8] and excitation transfer of biologi- 
cal system [9] also need to be treated as non-Markovian processes. Because of these distinct 
properties and extensive applications, more and more attention and interest have been de- 



voted to the study of non-Markovian processes o 
non- Markovianty 

aaor 



12 

23 



20 



and approaches 



21 



3JL the positivity 



14 



open systems, including the measures of 



16] , and some other dynamical properties 



22] of non-Markovian processes. Experimentally, the simulation 



24| of non-Markovian environment has been realized. 
The measure of non-Markovianty of quantum evolution is a fundamental problem which 
aims to detect whether a quantum process is non-Markovian and how much degrees it devi- 
ates from a Markovian one. So far, almost all measures can only serve to settle the former, 
i.e., serve as the sufficient (not necessary) condition for the emergence of non-Markovianty. 
Therefore, the problem for measuring the non-Markovianty of quantum processes still re- 
mains elusive and, in some sense, controversial. Based on the distinguishability of quantum 
states, Breuer, Laine and Piilo (BLP) [2] proposed a measure to detect the non-Markovianty 
of quantum processes which is linked to the flow of information between system and envi- 
ronment. Alternatively, Rivas, Huelga and Plenio (RHP) 10] also presented two measures 
of non-Markovianty by exploiting the dynamical behavior of quantum entanglement under 
the local trace-preserving CP maps. Other measures of non-Markovianty include the one 



proposed by Wolf et al. 
proposed by Lu et al. 



Ill ] based on the breakdown of the semigroup property, the one 
using quantum Fisher information flow, and the ones proposed 
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by Usha Devi et al. [12[ using relative entropy difference and fidelity difference. It is worth- 
while to stress that these measures of non-Markovianity are generally not the same. Thus, 
studying and exposing the relations among all these measures under some specific models 



becomes very important. Haikka et al. 



25| studied the links between the BLP measure 



and one of the RHP measures very recently, for a laser-driven qubit system embedded in a 
structured Lorentzian environment. By algebraic plus numerical calculations, they showed 
that, in the nonsecular regime, the two measures agree. But for other more cases, no definite 
result can be presented. 

In this paper, we study the equivalence of the three measures of non-Markovianity recently 
proposed by BLP and by RHP n]] in the case that a two-level system coupling to its 
environment via damped Jaynes-Cummings or dephasing model. Two important results are 
found: Firstly, the three measures are exactly equivalent in the aspect of detecting non- 
Markovianty for the involved models. And secondly, for our considered models, we find 
that the maximization in the definition of the BLP measure may be actually removed. This 
solves the problem that is being explored extensively by many researchers 26|, [271]. Our 
researches have the following features. First, the interaction models involved here are most 
fundamental in the theoretical studies of dynamics of open quantum systems, and we need 
not assuming any specific spectral density for the structured environment. Thus, our results 
possess good adaptability. Next, our operations are completely analytical and hence the 
results are more convincing. 

The article is organized as follows. In Sec. II we briefly review the three measures 
proposed by BLP and RHP respectively. The definition of equivalence between different 
measures is also presented. And in Sees. Ill and IV we study the equivalence of the three 
measures for a two-level system interacting with environment via damped Jaynes-Cummings 
model and dephasing model respectively. A conclusion is arranged in Sec.V. 



II. REVIEW OF THE THREE MEASURES 

In reference BLP presented a measure of non-Markovianty for quantum processes of 
open systems which is based on the idea that Markovian processes tend to continuously 
reduce the trace distance between any two states of a quantum system. Thus, an increase 
of trace distance during any time intervals implies the emergence of non-Markovianty. The 
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authors further linked the changes of trace distance to the flow of information between system 
and its environment, and concluded that the back flow of information from environment to 
the system is the key feature of a non-Markovian dynamics. Considering that the measure 
should reveal the total feature of a whole quantum process, they thus suggest the quantity, 



= max / dta(t, pi 2(0)), (1) 

Pl,2(0) J(T>0 

as the measure of non-Markovianty of a quantum process. Where cr(i, pi^O)) denotes the 
time derivative of trace distance for a pair of dynamical states with initial values pi,2(0) 
of the considered system. The time integration is extended over all intervals in which a 
is positive, and the maximum is taken over all pairs of initial states. For any Markovian 
process, M = 0. And if M > 0, the process must be non-Markovian. 

Alternatively, RHP [lO] presented another method to measure the non-Markovianty of a 
quantum process which is based on the monotonic drop of quantum entanglement between 
bipartite systems under the influence of local Markovian environments. Suppose a system 
of interest is initially prepared in a maximally entangled state with an ancillary particle, 
where only the system is influenced by a noise environment and the ancillary particle is free. 
Then the quantity for measuring the non-Markovianty of the quantum process is defined as 

X^ = f-\ dE[pS d f t)] \dt-^ (2) 

where the time derivative in the integrand is for the dynamical entanglement between the 
system and ancillary particle, AE = E[psA(to)] — E[psAit ma x)] denotes the difference of 
entanglements at the initial time to and the final time t max of the interest quantum process. 

The above two measures, Af and X^ E \ are introduced respectively through the mono- 
tonicity of trace distance or quantum entanglement under CP (or local CP) maps. Their 
mathematical forms are similar actually. It is not difficult to find that eq.(2) may be equiv- 
alently rewritten as 

X {E) = 2 f E (t)dt, (3) 

JE(t)>0 

with E(t) being the time derivative of dynamical entanglement E(t). Obviously, the measure 
is very similar to M in forms, with the time derivative E(t) of dynamical entanglement 
replaced by the time derivative a = D of the dynamical trace distance D(pi(t), p2{t)). The 
visible difference is: the measure M involves a maximum over all pairs of initial states, 
while I^ E > is defined via a given maximally entangled initial state and thus escapes from 



the optimization problem. Actually, the escaping of the optimization is just one of the 
original intentions for the measure to be proposed. In the following sections, we will 
demonstrate that, for the models considered in this article, the maximization process can 
be removed actually. 

The third quantity for measuring the non-Markovianty of a quantum process is also 
based on the entanglement dynamics between the system and an ancillary particle. Given 
the maximally entangled state |$) of the system plus the ancillary particle, a locally com- 
plete positive map e will keep the positivity of density operator p = ]$)($] invariable, i.e., 
> 0. Startin g fro m this point of view, a measure of non-Markovianty of quantum 



processes is then born [10] 



l = J o g(t)dt, (4) 

with 

= || [/ + (£, .8 J) M || - V 

v ' e->0+ e 

Where Ct is the super operator in the non-Markovian dynamical master equation % = C t (p) 
for the open system. 

The above three measures Af, and X, for the non-Markovianty of a quantum process 
are introduced through different ways. The physical meanings are different. A question 
naturally arises: They are equivalent each other? In this paper, we demonstrate that the 
three measures are equivalent each other when they apply to open two-level systems with 
damped Jaynes-Cummings or dephasing models. This equivalence implies that the measures, 
from different sides, well capture the intrinsical characters of non-Markovianty of quantum 
evolutional processes. 

As the end of this section, let us elaborate the implication for different measures to be 
equivalent. We say two measures Mi and M 2 to be equivalent means that: For a given 
quantum process, if by measure M±, the non-Markovianty emerges, then it also emerges by 
measure M 2 ; Conversely, if by measure M\, the non-Markovianty does not emerge, it also 
does not by measure M 2 . The exact logics are: Mi > -<=>- Af 2 > 0, Mi = -<=>- M 2 = 0. 
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III. THE CASE OF JAYNES-CUMMINGS MODEL 



In this section, we will demonstrate the equivalence of the three measures of non- 
Markovianty for the case of a two-level system coupled to its environment via Jaynes- 
Cummings model. For this end, we should first derive out the expressions of the three 
measures given the specific interaction model. For the calculation of TV, let us consider a 
two-level atom interacting with its environment via a Jaynes-Cummings model. The envi- 
ronment is assumed to be initially in the vacuum state, but its spectral density is arbitrary. 
This case can be solved exactly. For any pair of initial atomic states, pi(0) and /?2(0), one 
can obtain the time derivative of trace distance of the corresponding dynamical states asQ], 

, \ \ 2\G(t)\ 2 a 2 + \b\ 2 d 
f \G{t)\ 2 a 2 +\b\ 2dt 



where a = <1 |p x (0) 1 1> - <l|p 2 (0) [1> and b = <1 |p! (0) |0> - <l|p 2 (0)|0> are the differences of 
populations and of coherences respectively between the two given initial states. The function 
G(t) is defined as the solution of the integro-differential equation, 

^G{t) = - f dhf(t-h)G(t x ), (7) 



dt jo 

with initial condition G(0) = 1. The kernel f(t — ti) denotes the two-point reservoir corre- 
lation function which is the Fourier transformation of the spectral density. Introducing the 
time-dependent decay rate, 

70) = -2Re[||] = -^IGMI, (8) 
then the derivative of the trace distance may be rewritten as, 

*(*,Pi, 2 (0)) = - 7 (*)F(f). (9) 
Here the completely positive real function F(t) is defined as, 

a 2 e -|r(t) — I — ijj2 p -|r(t) 

F(t) = . + |b| (10) 
^a 2 e-nt) + [ftji 

with V(t) = $ dt'^(t'). According to eq.(l), the non-Markovian measure M hence may be 
expressed as, 

jV=-[ j(t)F(t)dt. (11) 

J^{t)<o 
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Where we assume that the pair of initial states ^1,2(0) just maximizes M so that the max- 
imization symbol is removed. Actually, due to the completely positivity of function F(t), 
the time intervals in which the trace distance increases monotonously, or equivalently the 
intervals in which non-Markovianty emerges (we call them non-Markovian intervals below), 
are uniquely determined by the condition j(t) < 0. The change of initial states pi,2(0) would 
not alter the positions and lengthes of these non-Markovianty intervals, i.e., would not alter 
the distribution of them. No matter how the initial states chnage, M would not shift from 
positive to zero, or vice versa. In this sense, the maximization to eq.(ll) may be removed 
without influencing the sensibility of M to detect non-Markovianty. This is an important 



result which solves the problem being explored by many researchers [26|, |27 ] . 

Next, we begin with the calculation of the second measure X^ E \ Consider a maximally 
entangled state of two two-level atoms, 

1$) * (|10) 5A +|01)<m), (12) 

where atom 1 is our considered system which couples to an environment via Jaynes- 
Cummings interaction and atom 2 is an ancillary particle which remains isolated from 
environment. The whole Hamiltonian reads, 



H = H Q + H I} (13) 

2 

i=l k 

Hi = ^gk^+ ] a k + h.c, 

k 

where and denote ladder operators for the i-th atom, Uk and are the frequency and 
annihilation operator for the /c-th harmonic oscillator of the environment, gu is the coupling 
constant. We assume the two atoms have the same transition frequency uq. Introducing to 
the interaction picture with respect Hq, one has, 

Hi = a\ 5 , fcafcexp[(w - u k )t] + h.c. (14) 

k 

Assume that the two atoms are initially in the entangled state |<&) and environment is 
initially in vacuum state |0), then the dynamical wave function for the compound system 
including both atoms and environment may be expressed as, 



m)) = Ci(t)\10) SA \0) E + c 2 (t)\01)sA\0} E + E Ck(t)\00} SA \l k } E . (15) 

k 

The dynamical reduced state for the two atoms thus reads, 

p SA (t) = (1 - | Cl (t)| 2 - |c 2 (t)| 2 )|00>(00| + | Cl (t)| 2 |10)(10| + |c 2 (t)| 2 |01>(01| 



+ci(t)^(t)|io)(oi| + c ;(*)p2(*)|oi)(io|, 

where the time-dependant coefficients are governed by Schrddinger equation, 



(16) 



-il\m)=Bi{t)W))- (17) 

After a straightforward deduction, we find c 2 (t) = c 2 (0) = 1/Vz and c\{t) = G(t)/\/2 
with G(t) being determined by eq.(7). If we use concurrence to describe the entanglement 
between the system and ancillary particle, we have from eq.(16) that, 

C(p S A(t)) = 2\cl(t)c 2 (t)\ = \G(t)l (18) 

and thus, 

jC{psA{t)) = -\l(t)e-^'\ (19) 
here the definition of ^(t) in eq.(8) is used. Eq.(5) then immediately produces, 

1^ = - [ 7 (t)e- r «/ 2 dt. (20) 

Jy(fi<0 

It again shows that the non-Markovian intervals for measure are uniquely determined 
by condition j(t) < 0. 

For the calculation of the third measure X, we need the dynamical master equation of 
the open two-level system lj], 

^ = ~S(t)[<T+(7^p]+ 7 (t)[<T-pa + - ±{<T+<T-, P }], (21) 

here S(t) = — 2Im(^y) with G(t) determined by eq.(7) and 7(t) is defined by eq.(8). After 
a straightforward calculation, one find from eq. (4) that 

. for 7 (t) > 
9(t) = { (22) 
-7(t) for 7(0 < 
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and thus 

T=-f i(t)dt. (23) 

Once again, the non-Markovian intervals for this measure are determined by 7(£) < 0. 

Through the above arguments, we conclude that, for an open two-level system with 
damped Jaynes-Cummings model, the distributions of non-Markovian intervals for the three 
measures, J\f, and X, are exactly the same, which are determined uniquely by the condi- 
tion 7(t) < 0. According to the viewpoint of equivalence for different measures mentioned in 
Sec. II, we thus conclude that the three measures are equivalent for detecting the existence of 
non-Markovianty of a two-level system interacting with environment via Jaynes-Cummings 
model. In the arguing process, we did not assume any specific structure of spectral density. 
Thus the result has good adaptability. They are valid for arbitrary structured environment 
(Lorentzian, Ohmic, waveguide spectrum etc.), arbitrary detunings (including resonance) 
between system and its environment, and arbitrary coupling strengthes. 



IV. THE CASE OF DEPHASING MODEL 



As the second example that we use to demonstrate the equivalence of measures of non- 
Markovianty, we consider a two-level atom which is coupled to a reservoir of harmonic 
oscillators via dephasing model. The Hamiltonian in Schrodinger picture is taken to be 

H = + E"* 6 * 6 * + E^( A 4 + Kh) (24) 
1 k k 

where ojq is transition frequency and a z — |1)(1| — |0)(0| the Pauli operator of the atom. u k , 
bk are respectively the frequency and annihilation operator for the fc-th harmonic oscillator 
of the reservoir. The coupling strength \ k is assumed to be complex in general. This 
dephasing model, which is extensively used to simulate the decoherence of a qubit coupling 
to its environment in quantum information science, can be solved exactly. For the initial 
state of the total system 

Psb(0) = p(0)®p B , (25) 

with p(0) the initial state of the atom and p B the initial thermal equilibrium state of the 
reservoir, the evolution for the elements of the reduced density matrix of the atom may be 
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written as []J 

Pu{t) = pn(0),p 00 (t) = poo(O) 



(26) 

Pio(t)=P* 01 (t) = p 10 (0)e r ^\ 
Here the negative dephasing function T p (t) is defined as 

roc 1 CD^djf 

TJt) = - / duJ{u)coth{u/2k B T) (27) 

JO OJ 

with J {to) the spectral density of the reservoir, ks the Boltzmann constant and T the 
reservoir temperature in thermal equilibrium. For any pair of initial states pi(0) and p 2 (0) 
of the atom, one can easily obtain the dynamical states p\{t) and p 2 (t), and then get the 
dynamical trace distance as, 



D{t, p lj2 (0)) = yja? + |6|2 e 2r p (t)_ (2 8) 

Where a = (l|p x (0)|l) - (l|p 2 (0)|l) and b = (l\ Pl {0)\0) - (l|p 2 (0)|0) are the differences of 
the populations and of coherences for the two given initial states respectively. Inserting this 
trace distant into eq.(l), one immediately get the non-Markovian measure for this dephasing 
model as 

r 1512 2r p (i) 

M p = -2 dt lp {t) r ±L- (29) 

J Mt) <o ^J a 2 + |6|2 e 2r p (t) 

where the dephasing rate is defined as j p {t) = — ~T p {t), and based on the same reason as 
before, we neglect the maximization to M p . It shows that the non-Markovian intervals are 
uniquely determined by j p {t) < 0. 

In order to calculate the expression of the second measure we introduce, in like 

manner as before, an isolated ancillary atom which is initially prepared in a maximally 
entangled state |$) in eq.(12) with the system atom. According to the evolution of eq.(26) 



and employing the trick in reference 
for the system and ancillary atoms, 



one can easily obtain the dynamical density matrix 



p 5A (t) = ^[|10)(10|+e r ^|01)(10| + e r ^|10)(01| + |01)(01|]. (30) 

Where the first state refers to the system and the second one to the ancillary atom. It is 
easy to find that the concurrence for this state is C(t) = exp[T p (t)], and thus eq.(5) produces 

Xf ) = -4 / lp (t)e r ^dt. (31) 
h P (t)<o 
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Again, the non-Markovian intervals for measure X^ is also uniquely determined by j p (t) < 
0. 

Finally, let us calculate the measure X of eq.(3). It is not difficult from eq.(26) to verify 
that the dynamical master equation for the dephasing two-level atom can be written as 

-£ = 7 P (t)[a z p(T z -p], (32) 

where we again have used the relation T p (t) = —2j p (t). Following the same deduction for 
eqs.(21)-(23), we obtain 

X p = -2 / lp {t)dt. (33) 

Once again, the non-Markovian intervals for this measure are uniquely determined by j p (t) < 
0. 

All in all, for a two-level system with dephasing model, we again draw the conclusion that 
the distributions of non-Markovian intervals for the three measures are exactly the same and 
uniquely determined by the condition j p (t) < 0. Thus the three measures are equivalent. 



V. CONCLUSION 



In conclusion, we have compared the three measures of non-Markovianty proposed in 
references 3] and [h]] respectively for two-level systems interacting with environments via 
damped Jaynes-Cummings or dephasing model. The results show that the three measures 
have exactly the same distribution of non-Markovian intervals, and hence are actually equiv- 
alent each other in the aspect of detecting non-Markovianty of quantum processes. This 
equivalence implies that these measures, in different ways, well capture the intrinsical char- 
acters of non-Markovianty of quantum processes. We have also found that the maximization 
in the BLP measure can be removed for the considered models, without influencing the sen- 
sibility of the measure to detect non-Markovianty. This result is important which extremely 
avoids the complicated mathematical calculations. The dynamical models considered here, 
i.e., damped Jaynes-Cummings and dephasing models, represent two kinds of fundamental 
coupling forms in studying problems of open systems. And we did not assume any specific 
spectral density for the structured environment. Thus our results have good adaptability. 
They apply to arbitrary structured environment, arbitrary detunings between system and 
its environment, and arbitrary coupling strengthes. Exposing the relations between different 

11 



measures of non-Markovianty in different physical systems and different models is impor- 
tant. We expect similar researches for many more measures could be expanded for more 
general models or in more complicated systems. 
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